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1. Introduction 



Strong coupling and quark masses constitute fundamental parameters of QCD part in the stan- 
dard model. One of the important task of lattice QCD is to determine these parameters from inputs 
at low energy scale. Hadron masses, meson decay constants and quark potential quantities are 
adopted as physical inputs and QCD running coupling and quark masses should be given. These 
results may be compared with independent evaluation from high energy inputs, which may uncover 
existence of systematic error if there is. 

In a course of evaluating these fundamental parameters in the Lagrangian we need a process of 
renormalization in some scheme. It is now recognized that systematic deviation due to perturbative 
renormalization is large for quark mass and also for running coupling at low energy region. Non- 
perturbative renormalization is essential for this work. Among several non-perturbative schemes 
on the lattice the Schrodinger functional (SF) scheme ||, ||, has an advantage that systematic 
errors can be unambiguously controlled: A unique renormalization scale is introduced through 
the box size to reduce the lattice artifact and a large range of the renormalization scale can be 
covered by the step scaling function (SSF) technique. The latter virtue matches our purpose to 
make comparison with high energy inputs. 

For the SF scheme we need to start with evaluation of the running coupling in order to in- 
troduce the renormalization scale. It is strongly expected that three light quarks should involve 
in running of the coupling at low energy scale, where non-perturbative effects becomes large. So 
introduction of u,d, s quarks would be important for non-perturbative running of the coupling. This 
report presents preliminary results of our calculation for the running coupling in Nf = 2 + 1 QCD 
with SF scheme. The physical scale shall be introduced through the Sommer scale ro evaluated 
independently by the CP-PACS collaboration with light Nf = 2 + 1 configuration [|5p. 



2. Schrodinger functional formalism and action 

The Schrodinger functional is given as a field theory in a finite box of size L 4 with a Dirichlet 
boundary condition at temporal boundary. For QCD the Dirichlet boundary condition is set for 
spatial component of the gauge link 

• , \ 

U k (x) \ X0=0 = exp (aC k ) , U k (x) \ Xo=T = exp (aC' k ) , C { k ] = ~ <j>P (2. 1) 

and quark fields 

W(x) |xo=0 = V(*) \x =T = 0, ljf(x) |xo=0 = WW \x =T = 0. (2.2) 
Under a lenient condition it is proved that the tree level gauge effective action has a global minimum 



around a background field which is uniquely given by the boundary fields (2.1). On the other 
hand the fermionic mode is shown to have a mass gap, with which we are able to define mass 
independent scheme in the chiral limit without any extrapolation. The renormalization scale is 
given only by the box size L. 
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We adopt the renormalization group improved gauge action of Iwasaki type 

S s = l I W (C,^)Retr(l-P(C)) + £ £ m(C,^)Retr(l-/?(C)), (2.3) 

iv CeS v CeSi 

where So and 5i are sets of oriented plaquettes and rectangles. The weight factor Wq/\ is given to 
cancel the 0(a) contribution from the boundary according to @, 0]. The boundary improvement 
coefficients are set to tree level value cf = 1 and cf = 3/2, which is shown to give better scaling 
behavior than one loop value We take the same values for boundary link (2T) as in the 

previous work of the Alpha collaboration [||, §]. 

We used the improved Wilson fermion action with clover term 

Sf = a 41 £W(Dw + m Q )Y, D w = - (ty (V M + V*) -aV* V M ) -csw-ff^y. (2.4) 
x Z • 4 

The improvement coefficient csw is given non-perturbatively in a polynomial form [^] which covers 
1-9 < j8 < 12.0. We notice that there is a contribution from the boundary to cancel the 0(a) effect 
there 

S (a) = a 3 ^(^-l)(W(x,l)Y(x,l)+W(x,T-y)w(x,T-y)), (2.5) 

x 

for which the one loop value JTQ| ] is taken c f = 1 — 0.0088 l(28)gg. We set twisted periodic bound- 
ary condition in spatial direction y(x + Lk) = e l6 \j/( x ) with 9 = n/5 @ §]. 

The renormalized gauge coupling in the SF scheme is defined as a coefficient of the effective 
action r[V^] at the global minimum. For numerical simulation we take derivative in terms of a 
parameter r\ introduced in the background field 0, and define the SF coupling as pi] 



i dr[v^ 



t(L) k dri 
where k is a normalization coefficient evaluated at tree level 



(2.6) 

77=0 



3. Our strategy 

The goal of our project for the running coupling is to derive the renormalization group invariant 
(RGI) scale Aqcd in a unit of the Sommer scale r^. The RGI scale A is scheme dependent and is 
defined as follows for the SF scheme 

Asp = I (WMr*«p (-^y exp (- [\ s (^ + _L,_ £)) , (3,) 

where g(L) is a renormalized coupling in SF scheme at a scale L and j3 (g) is renormalization group 
function (/3 -function) with its perturbative expansion coefficients 

JS (g) = ~g 3 (bo + big 2 + b 2 g 4 + •••)• (3.2) 
Derivation of the RGI scale is given by the following steps in the SF scheme. 
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(i) We start by calculating the SSF £ on the lattice at several box sizes and lattice spacings. The 
SSF gives a relation between the renormalized couplings when the renormalization scale is changed 
by factor two Y,(u,a/L) = g 2 (2L)\ u __ 2 ^ [§01, where the scale is given by the renormalized cou- 
pling g 2 (L) and discretization error by a/L. Taking the continuum limit o(u) = lim a / L ^ L(u,a/L) 
and performing a polynomial fit we have a full non-perturbative running of the coupling in a dis- 
cretized manner. 

(ii) In the second step we define a reference scale L max through a fixed value of renormalized 
coupling g 2 (L max ). The value of g 2 (L max ) is rather ambiguous if it is well in low energy region. We 
then start from L max and follow non-perturbative RG flow through the SSF into high energy region. 
After n ~ 8 iterations the scale L = 2~"L max is already in perturbative region where discrepancy 
between perturbative and non-perturbative RG running is negligible. 

(iii) Substituting g 2 (L) and L = 2~ n L max given in the above into (3T) and evaluating the inte- 
gral with three loops /3 -function in the SF scheme [|l|] we get the RGI scale AsF^max in terms of 
the reference scale. 

(iv) In the last step we need the physical input ro measured in an independent large scale 
simulation at some lattice spacing a. The reference scale should also be measured at the same 
lattice spacing to give a ratio ro/L max . The requirement for the lattice spacing and the reference 
scale is that magnitude of the lattice artifact a/ro and a/L max should be kept small. Multiplying 
these factors we get the RGI scale Asf?"o in terms of the Sommer scale. Transformation into the 
MS scheme is given exactly at one loop A^ = 2.612Asf- 

4. Step scaling function 

We adopted seven renormalized couplings to cover from the weak coupling region g 2 = 1.001 
to strong region g 2 = 3.418 separated approximately by twice the renormalization scale. For each 
coupling we used three boxes L/a = 4, 6, 8 to take the continuum limit. 

HMC algorithm is adopted for two flavours and RHMC algorithm for the third flavour, all of 
which are taken to be the common mass. We adopted CPS++ code and modified for SF formalism. 
For machines we make use of T2K, PACS-CS and PC cluster Kaede at University of Tsukuba, T2k 
and SRI 1000 at University of Tokyo and PC cluster RSCC at Riken. 

We start by tuning the value of j8 and K to reproduce the same renormalized coupling at 
each box sizes keeping the PCAC mass to zero, where the PCAC relation is defined in terms of 



the improved axial current with non-perturbative improvement coefficient [12]. We notice that 
distribution of inverse of the coupling 1 /g 2 turned out to be smooth Gaussian even at the lowest 
energy scale [JsJ] . This is contrary to the standard Wilson gauge action [|J] and we need no re- 
weighting. 

The renormalized coupling g 2 (2L) at larger scale is modified perturbatively in order to cancel 
deviation of the PCAC mass from zero and that of the renormalized coupling g 2 (L) from a fixed 
value[|, U, 01. A part of 0{ a) error is canceled at one loop level with coefficients given in 
Ref. In the end we get the 0(a) improved SSF on the lattice. 

Preliminary result is plotted in figure |TJ The left panel shows scaling behavior of the SSF 
at each renormalization scale, which turned out to be good except at the strongest coupling. We 
performed three types of continuum extrapolation: constant extrapolation with finest two (filled 
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Figure 1: SSF on the lattice with its continuum extrapolation at each renormalization scale (left). RG flow 
oftheSSF (right). 



-0.04 




Figure 2: Non-perturbative j3 -function for Nf = 3 and 2 QCD. 



symbols) and three data points (open symbols) and linear extrapolation with three data (open cir- 
cles). As is plotted in the figure they are consistent with each other except at the strongest. We 
adopted the continuum fit with finest two lattice spacings as our preliminary continuum value. In 
the right panel the RG running of the SSF is plotted. We divide the SSF with the coupling g 2 (L) to 
get better resolution. Polynomial fit of the continuum SSF to sixth order 

a(ll) = U + SqU 1 + SiM 3 +S2U 4 + Sjll 5 + S4U, 6 (4.1) 

is plotted (solid line) together with the three loop perturbative running (dotted line). We used one 
and two loop values for so and *i in the fit. From the polynomial form of the SSF we derive the 
non-perturbative /3-function of the Nf = 3 QCD, which is plotted in figure ^[ The /3-function of 
Nf = 2 QCD is reproduced from data of the Alpha collaboration [Q] for comparison. 

5. Introduction of physical scale 

From a independent simulation of PACS-CS collaboration at j3 = 1.90 ^ preliminary value 



5 



Non-perturbative renormalization ofNf = 2+1 QCD with SF scheme 



Yusuke Taniguchi 



of the Sommer scale is given as a/ro = 0.131(35) in the chiral limit. Evaluation of the strong 
coupling in the SF scheme at the same j3 in 4 4 box gives g 2 (L max ) = 4.695(23) in the chiral limit. 
We adopt this coupling as a definition of L max . 1 

Starting from w max = 4.695(23) we iterate non-perturbative RG flow eight times according to 
polynomial fit (14) and substitute the result into (34). In the end we get Asr^max = 0.238(19). 
Multiplying a/ro and L max /a = 4 we change the reference scale to r$ and we have Asf?o = 
0.45(12). The large error mainly comes form systematic error of a/ro during chiral extrapola- 
tion of strange quark mass. We need few more points to take the massless limit in a rigid way. In 
order to give Aj^§- in a unit of MeV we also need to check validity of ro = 0.5 fm in the chiral limit. 



6. Step scaling function for quark mass 

Since we are calculating the PCAC mass with the 0(a) improved axial current it is possible 
to derive renormalization factor for the pseudo scalar density as a byproduct, from which we can 
extract non-perturbative running of the quark mass. However the scaling behavior of the pseudo 



scalar density was shown to be bad even perturbatively [14] under the inhomogeneous boundary 



gauge field (2.1) and twist factor 9 = n/5 in spatial direction. 



In this report we define the pseudo scalar density renormalization factor as 

Zp(go,L/a) = — — — — , (6.1) 



'(lattice) 

where propagators fp and f\ are given in [gj. We expect cancellation of 0(a) effect at tree level 
dividing by tree level propagator on the lattice. The SSF on the lattice is given by 



a\ _ Zp(g ,2L/a) 

Li . 



L) Zp(g Q ,L/a) 



(6.2) 

g 2 (L)=u,m=0 



The result is plotted in figure [| which shows a rather bad scaling behavior. Although we take the 
continuum limit with finest two lattice spacings in this report, we may need finer lattice spacings 
to reduce systematic uncertainty, which seems to be unrealistic for computational cost. 



7. Conclusion 

We present a preliminary result for the Nf = 2+1 QCD running coupling in the mass indepen- 
dent SF scheme in the chiral limit. We used seven scales to cover from low energy to high energy 
region and three lattice spacings to take the continuum limit at each scale. 

Tuning of j8 and K has been completed to fix seven scales in the massless limit. We are 
now evaluating the SSF at the finest lattice spacings. Our preliminary result shows a good scaling 
behavior except at the lowest energy scale, for which we may need finer lattice spacing to take the 
continuum limit. In order to evaluate Aj^ precisely we need to derive a/ro in the chiral limit in 
more rigid way. Main source of the systematic error is an extrapolation of the strange quark mass 

'Unfortunately simulation in 6 4 box gives g 2 (L) = 6.71(16), which exceeds our largest coupling g 2 = 5.35(10) for 
the SSF. 
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Figure 3: SSF of the pseudo scalar density on the lattice with its continuum extrapolation at each renormal- 
ization scale (left). RG flow of the SSF (right). 



and we are planning to perform simulation at different parameters for r$. We also need to check 
validity of r$ = 0.5 fm in the chiral limit before we evaluate A^- in terms of MeV. 

The scaling of the pseudo scalar density SSF is rather bad under inhomogeneous background 
gauge field. We may need better setup with vanishing background field, which we are planning as 
a next step. 

This work is supported in part by Grants-in-Aid of the Ministry of Education, Culture, Sports, 
Science and Technology-Japan (NOs. 18740130, 18104005). 
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